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ABSTRACT 

Aims.  We  present  a  hot  white  dwarf  (WD)  luminosity  function  (LF)  using  data  taken  from  the  Sloan  Digital  Sky  Survey  (SDSS) 
Data  Release  4.  We  present  and  discuss  a  combined  LF,  along  with  separate  DA  and  non-DA  as  LFs.  We  explore  the  completeness 
of  our  LFs  and  interpret  a  sudden  drop  in  the  non-DA  LF  near  2  Mboi  as  a  transition  of  the  non-DA  WD  atmosphere  into  the  DA  one 
during  WD  evolution.  Our  LF  extends  roughly  between  -0.5  <  Mboi  <  7  or  equivalently,  ~  120  000  K  >  Teg  >  ~25  000  K.  Our  LF 
should  now  be  useful  for  estimates  of  recent  star  formation  and  for  studies  of  neutrino  and  other  potential  particle  emission  losses  in 
hot  WDs. 

Methods.  To  create  a  sample  whose  completeness  can  be  characterized  fully,  we  used  stars  whose  spectra  were  obtained  via  the 
SDSS’s  “hot  standard”  target  selection  criteria.  The  hot  standard  stars  were  purposefully  targeted  to  a  high  level  of  completeness  by 
the  SDSS  for  calibration  purposes.  We  are  fortunate  that  many  of  them  are  hot  white  dwarfs  stars.  We  further  limited  the  sample  to 
stars  with  fitted  temperatures  exceeding  23  500  K  and  log  g  >  7.0.  We  determined  stellar  distances  for  our  sample  based  on  their 
absolute  SDSS  g  filter  magnitudes,  derived  from  WD  stellar  atmosphere  model  fits  to  the  SDSS  stellar  spectra. 

Results.  We  compared  our  LF  with  those  of  other  researchers  where  overlap  occurs;  however,  our  LFs  are  unique  in  their  extension 
to  the  most  luminous/hottest  WDs.  The  cool  end  of  our  LF  connects  with  the  hot  end  of  previously  determined  SDSS  WD  LFs  and 
agreement  here  is  quite  good.  It  is  also  good  with  previous  non-SDSS  WD  LFs.  We  note  distinct  differences  between  the  DA  and 
non-DA  LFs  and  discuss  the  reliability  of  the  DA  LF  at  its  hot  end.  We  have  extended  the  range  of  luminosities  covered  in  the  most 
recent  WD  LFs.  The  SDSS  sample  is  understood  quite  well  and  its  exploration  should  contribute  to  a  number  of  new  insights  into 
early  white  dwarf  evolution. 

Key  words,  stars:  luminosity  function,  mass  function  -  white  dwarfs 


1.  Introduction 

The  luminosity  function  (LF)  and  space  density  of  WDs  pro¬ 
vide  constraints  on  the  local  star  formation  rate  and  history 
of  the  Galactic  disk  in  the  Solar  neighborhood.  They  can  also 
be  used  to  help  determine  the  role  of  weakly  interacting  parti¬ 
cles  in  WD  cooling  (Irsen  et  al.  2008)  and  thus  their  ages.  A 
large  WD  sample  of  known  completeness  is  required  to  deter¬ 
mine  a  precise  LF  to  address  these  questions.  Luckily,  the  Sloan 
Digital  Sky  Survey  (SDSS)  (Gunn  et  al.  1998,  2006;  York  et  al. 
2000;  Stoughton  et  al.  2002;  Pier  2003;  Tucker  et  al.  2006)  pro¬ 
vides  the  perfect  sample  for  this  effort.  Although  primarily  de¬ 
signed  as  an  extra-Galactic  survey,  the  SDSS  is  an  extremely 
rich  source  of  stellar  (and  in  particular,  white  dwarf)  spectra. 
Eisenstein  et  al.  (2006a)  published  a  catalog  of  spectroscopic¬ 
ally  confirmed  white  dwarf  stars  from  the  SDSS  Data  Release  4 
(DR4:  Adelman-McCarthy  et  al.  2006).  Their  catalog  contained 
9316  WDs,  nearly  3.6  times  the  number  of  WDs  found  in  the 


first  SDSS  WD  catalog  (Kleinman  et  al.  2004)  based  on  DR1 
(Abazajian  et  al.  2003)  data,  consistent  with  the  ratio  of  spec¬ 
troscopic  sky  coverage  (4783  /  1360  square  degrees)  in  DR4 
vs.  DR1. 

Harris  et  al.  (2006)  have  constructed  a  reduced  proper- 
motion  based  WD  LF  using  the  photometric  sample  from 
SDSS  DR3  (Abazajian  et  al.  2005)  and  proper  motions  from 
Monet  et  al.  (2003).  The  large  data  set  available  from  SDSS 
enabled  them  to  sample  the  LF  within  a  magnitude  range  of 
~7  <  Mboi  <  14  using  0.5  mag  bins,  with  each  bin  containing 
hundreds  of  stars  for  quite  excellent  statistics.  Hu  et  al.  (2007) 
also  calculated  a  DA  WD  LF  based  on  the  much  smaller  SDSS 
DR1  WD  sample,  but  extending  to  Mboi  =  4  at  its  blue  end. 
These  LFs  compare  well  with  that  of  Liebert  et  al.  (2005),  which 
is  based  on  DAs  from  the  Palomar  Green  survey  (Green  et  al. 
1986).  Most  recently,  DeGennaro  et  al.  (2008)  have  presented 
a  DA  WD  LF  also  based  on  SDSS  DR4  data,  using  statisti¬ 
cal  arguments  and  a  comparison  to  the  Harris  et  al.  (2006)  LF 


Article  published  by  EDP  Sciences 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2009  to  00-00-2009 


4.  TITLE  AND  SUBTITLE 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


A  hot  white  dwarf  luminosity  function  from  the  Sloan  Digital  Sky  Survey 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Mt.  Suhora  Observatory, Cracow  Pedagogical  University,  report  number 

ul,Podchorazych  2,  30-084  Cracow,  Poland, , 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Aims.  We  present  a  hot  white  dwarf  (WD)  luminosity  function  (LF)  using  data  taken  from  the  Sloan 
Digital  Sky  Survey  (SDSS)  Data  Release  4.  We  present  and  discuss  a  combined  LF,  along  with  separate  DA 
and  non-DA  as  LFs.  We  explore  the  completeness  of  our  LFs  and  interpret  a  sudden  drop  in  the  non-DA 
LF  near  2  Mbol  as  a  transition  of  the  non-DA  WD  atmosphere  into  the  DA  one  during  WD  evolution.  Our 
LF  extends  roughly  between  &#8722;0.5  <  Mbol  <  7  or  equivalently,  &#8764;120  000  K  >  Teff  > 

&#8764;25  000  K.  Our  LF  should  now  be  useful  for  estimates  of  recent  star  formation  and  for  studies  of 
neutrino  and  other  potential  particle  emission  losses  in  hot  WDs.  Methods.  To  create  a  sample  whose 
completeness  can  be  characterized  fully,  we  used  stars  whose  spectra  were  obtained  via  the  SDSS?s  ?hot 
standard?  target  selection  criteria.  The  hot  standard  stars  were  purposefully  targeted  to  a  high  level  of 
completeness  by  the  SDSS  for  calibration  purposes.  We  are  fortunate  that  many  of  them  are  hot  white 
dwarfs  stars.  We  further  limited  the  sample  to  stars  with  fitted  temperatures  exceeding  23  500  K  and  logg 
>  7.0.  We  determined  stellar  distances  for  our  sample  based  on  their  absolute  SDSS  g  filter  magnitudes, 
derived  from  WD  stellar  atmosphere  model  fits  to  the  SDSS  stellar  spectra.  Results.  We  compared  our  LF 
with  those  of  other  researchers  where  overlap  occurs;  however,  our  LFs  are  unique  in  their  extension  to 
the  most  luminous/hottest  WDs.  The  cool  end  of  our  LF  connects  with  the  hot  end  of  previously  determined 
SDSS  WD  LFs  and  agreement  here  is  quite  good.  It  is  also  good  with  previous  non-SDSS  WD  LFs.  We  note 
distinct  differences  between  the  DA  and  non-DA  LFs  and  discuss  the  reliability  of  the  DA  LF  at  its  hot  end. 
We  have  extended  the  range  of  luminosities  covered  in  the  most  recent  WD  LFs.  The  SDSS  sample  is 
understood  quite  well  and  its  exploration  should  contribute  to  a  number  of  new  insights  into  early  white 
dwarf  evolution. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


340 


J.  Krzesinski  et  al.:  Hot  white  dwarf  luminosity  function 


to  estimate  sample  completeness.  The  agreements  and  disagree¬ 
ments  between  the  two  LFs  are  fairly  well  understood,  although 
their  cutoff  at  Mboi  ~  6  offers  only  a  small  chance  of  compar- 
ision  to  our  own  LF,  which  extends  to  much  hotter  stars.  Prior 
to  the  flurry  of  new  SDSS-based  WD  LFs,  WD  LFs  were  con¬ 
structed  by  for  example  Flemming  et  al.  (1986)  and  Vennes  et  al. 
(2002).  There  were  also  attempts  to  estimate  the  local  non-DA 
WDs  space  density  dating  back  to  1985  (Wesemael  et  al.  1985). 
None,  however,  extend  below  Mb0i  =  4  mag,  while  our  new  LF 
extends  to  below  Mboi  =  -0.5  mag. 

Traditionally,  distances  in  WD  LFs  have  been  calculated  us¬ 
ing  a  statistical  proper  motion  method  (e.g.  Liebert  et  al.  1988; 
Evans  1992;  Oswalt  et  al.  1996,  and  others).  For  our  relatively 
bright,  therefore  distant  stars,  this  approach  is  less  likely  to  work, 
since  many  are  simply  too  far  away  to  have  measurable  proper 
motions.  We  therefore  calculated  bolometric  magnitudes  derived 
from  fitting  theoretical  atmosphere  model  to  the  SDSS  spectra 
and  compared  them  to  the  observed  SDSS  g  filter  (SDSS  filters: 
Fukugita  et  al.  1996)  magnitudes  to  obtain  distances  to  each  of 
our  objects.  Unlike  all  the  other  WDLFs,  ours  is  based  solely  on 
spectroscopically-confirmed  WDs  from  the  SDSS  spectroscopic 
sample  and  covers  a  Mboi  range  from  -0.75-7.0. 

2.  Our  hot  white  dwarf  sample 

The  SDSS  “hot_std”  category  targets  bright  (g  <  19),  isolated, 
hot  stars  for  calibration  purposes.  In  short,  the  following  crite¬ 
ria  define  the  hot_stds,  where  the  subscript  o  indicates  a  fully- 
dereddened  magnitude: 

-  g  >  14 

-  ga  <  19 

-  -1.5  <  (m  -  g)0  <  0 

- 1.5  <  (g  -  r) o  <  0. 

Because  of  the  relative  rarity  of  these  objects  in  the  sky  and  their 
importance  in  calibration,  the  SDSS  gave  hot  standard  stars  the 
highest  priority  in  assigning  fibers  (Blanton  et  al.  2003).  In  the¬ 
ory,  then,  WDs  within  the  color  selection  range  of  the  hot_std 
targets  could  be  100%  complete  within  the  SDSS  database. 
However,  because  the  hot_stds  are  used  for  calibration  some 
fairly  severe  additional  non-photometric  requirements  are  placed 
on  them  in  order  that  the  photometry  and  spectroscopy  can  be 
completely  trusted.  These  additional  constraints  require  the  tar¬ 
gets  to  have  no  interpolated  pixels  and  to  be  unblended  and  far 
from  CCD  edges.  Eisenstein  et  al.  (2006a)  estimated  that  these 
additional  constraints  lower  the  spectroscopic  completeness  of 
the  photometrically-specified  sample  to  below  77%,  contrary  to 
the  earlier  statements  of  Krzesinski  et  al.  (2004).  Nonetheless, 
the  hot_stds  are  still  among  the  most  complete  of  all  stars  spec¬ 
troscopically  observed  by  the  SDSS  and  their  incompleteness 
is  calculable  (see  Sect.  3.2).  Since  the  Eisenstein  et  al.  (2006a) 
WD  catalog  does  not  state  which  of  their  WDs  were  targeted 
as  hot_stds,  we  supplemented  their  catalog  with  SDSS  DR4 
database  queries  for  all  objects  targeted  as  hot_stds. 

Figure  1  is  a  plot  of  hot  DAs  and  DBs  from  Eisenstein  et  al. 
(2006a)  along  with  the  hot_std  (u  -  g)0  and  (g  -  r)a  color  cuts. 
The  (m  -  g)Q  <  0.0  color  selection  ends  up  cutting  out  DAs  at 
cooler  temperatures  than  the  DBs.  Objects  above  the  horizontal 
line  (some  DAs  cooler  than  23  500  K)  are  not  targeted  as  hot_std 
stars,  thereby  introducing  some  incompleteness  in  our  sample 
below  that  temperature.  We  therefore  imposed  an  additional  cut 
on  our  sample  of  77 fr  >  23  500  (the  vertical  line  in  Fig.  1),  so 
the  temperature  selection  of  DAs  and  non-DAs  in  our  sample  is 
unaffected  by  the  DAs  higher  ( u  -  g)  values. 


Fig.  1.  The  photometric  ( g  -  r)„  and  (u  -  g)a  (small  and  large  sym¬ 
bols  respectively)  vs.  fit  temperatures  for  our  hot_std  sub-sample  of  the 
DR4  WD  stars  from  Eisenstein  et  al.  (2006a).  Only  DA  (solid  squares) 
and  DB  (open  circles)  WDs  with  Tc g  >  20  000  K  are  shown.  The  dashed 
line  shows  our  additional  TeS  sample  cut  at  23  500  K. 


3.  Luminosity  function 

Once  our  white  dwarf  sample  was  well-defined,  we  determined 
their  distances  and  absolute  bolometric  magnitudes.  Next,  we 
calculated  the  sample’s  completeness  and  volume  corrections 
before  forming  our  final  luminosity  functions. 

3.1.  The  atmospheric  models  and  distances 

We  used  the  Eisenstein  et  al.  (2006a)  fit  Teg  and  log  g  parame¬ 
ters  for  the  pure  DA  and  DB  stars  in  our  sample.  Unfortunately, 
Eisenstein  et  al.  (2006a)  do  not  include  any  predicted  absolute 
photometric  magnitudes  in  their  output,  so  we  needed  additional 
work  to  convert  the  77ff  and  log  g  fits  into  an  absolute  magni¬ 
tude  which  we  could  use  to  calculate  distances.  Bergeron  (2001, 
private  communication)  calculated  a  grid  of  log  g,  7'cfr  vs.  Mbol 
and  absolute  SDSS  u,  g,  r.  i,  z.  filters  magnitudes  for  pure  DA 
and  DB  WDs.  Eisenstein  et  al.  (2006a)  demonstrated  that  their 
models  (from  Detlev  Koester;  published  in  Finley  et  al.  1997) 
are  in  good  agreement  with  Bergeron’s,  so  we  can  use  the  Teg 
and  log  g  fits  from  their  Koester  model  fits  to  obtain  absolute 
SDSS  filter  magnitudes  from  the  Bergeron  grid.  We  used  these 
absolute  magnitudes  to  compare  with  the  observed  magnitudes 
and  obtain  distances  to  all  our  objects. 

For  the  non  DA/DB  stars,  we  used  other  spectral  models 
and  fits  accordingly:  Rauch  &  Deetjen  (2003)  non-LTE  H+He 
models  for  the  DAOs  and  Hiigelmeyer  et  al.  (2006)  models 
for  DO  and  PG1159  spectra.  We  used  stellar  radii  taken  from 
evolutionary  models  of  Wood  (1995)  with  thick  and  thin  hy¬ 
drogen  layers  for  DAO  and  DO  WDs  and  Bloecker  (1995), 
Schoenberner  (1983)  and  Wood  &  Faulkner  (1986)  models  for 
PG1 159  stars.  DAB  (3  stars)  and  DBA  (1  star)  WDs  were  treated 
in  the  same  way  as  DA  and  DB  WDs,  respectively.  For  magnetic 
or  binary  DAs  and  DBs  the  fit  parameters  were  taken  from  the 
Eisenstein  et  al.  (2006a)  DR4  WD  catalog;  for  binary  DOs,  from 
Hiigelmeyer  et  al.  (2006).  Finally  we  calculated  two  LFs,  one  for 
all  WDs  including  magnetic  and  binary  WDs,  and  one  without 
them.  The  LFs  with  and  without  the  magnetic  and  binary  WDs 
are  not  significantly  different  from  each  other,  so  we  therefore 
continue  with  these  objects  included  in  our  further  analysis. 

In  general,  an  absolute  SDSS  filter  magnitude  A7n,  can  be 
obtained  by  implementing  the  modified  Girardi  et  al.  (2002) 
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formula  (Eq.  (2)): 

Mflt=  2.5  log  (3631  x  1(T23)  +  2.5  log  ^  '  J  c  { dA 
+  5  log  pc  •  10  -  5  log(7?0  ■  R*) 

-  2.5  log  |;r  ■  A  ■  A2  cm  ■  Flux  ■  Filter! A)  dAj 

+  zeropoint(flt) 
where: 


fit 

3631  x  1(T23  [erg/cm2/s/cm] 


c 

A 2  cm 

pc  =  3.0857  x  1018  [cm] 
Re 

Flux  ■  n  [erg/cm2/s/cm] 
Filter(A) 

A  [A] 

R+  [Re] 

zeropoint(flt)  [mag] 


one  of  u,  g.  r,  i,  z  SDSS  filters, 
AB  system  zero  magnitude  for 
a  SDSS  filter 
Oke  &  Gunn  (1983), 
speed  of  light  in  cm/s, 
Angstrom  to  cm  conversion, 

1  parsec, 

solar  radius  in  cm, 
flux  density  from  the  model, 
SDSS  filter  transmission  within 
the  wavelength  range  A1-A2, 
wavelength  in  Angstroms, 
stellar  radius  taken  from 
evolutionary  models 
SDSS  photometric  zero  point. 


Here  zeropoint(u)  =  -0.04,  zeropoint(gr)  =  zeropoint(r)  =  0.00, 
zeropoint(i)  =  0.015  and  zeropoint(z)  =  0.03.  The  SDSS  magni¬ 
tude  system  (photometric  system:  Smith  et  al.  2002;  Ivezic  et  al. 
2004;  Hogg  2001)  is  close  to  being  on  the  AB  system  Oke  & 
Gunn  (1983),  but  not  quite,  so  the  above  zero  points  are  needed 
(Abazajian  et  al.  2004;  Eisenstein  et  al.  2006b). 

We  used  the  Girardi  et  al.  (2002)  equations,  above,  to  obtain 
absolute  SDSS  g  filter  magnitudes  for  all  WDs  in  our  sample. 
For  most  of  them,  we  calculated  absolute  bolometric  magnitudes 
(by  extending  the  wavelength  range  in  the  equation  from  0  to  in¬ 
finity).  For  some  of  the  DOs  and  PG1159  stars  in  the  sample, 
however,  we  do  not  have  complete  enough  models  to  estimate  a 
total  stellar  flux,  so  in  those  cases  we  used  the  standard  equation, 
below,  to  obtain  absolute  bolometric  magnitudes: 

Mbo,  =  42.39  -  5  ■  log  /?*  -  10  ■  log  reff . 

We  compared  results  using  this  equation  with  the  absolute  mag¬ 
nitudes  determined  by  our  model  grids  and  found  the  answers  to 
be  in  good  agreement  (~0.3  mag  on  average). 

We  do  not  have  models  for  the  hottest  DBs  in  our  sample 
those  between  about  30000  K  and  42  500  K,  the  so-called  “DB 
gap”  region  (Eisenstein  et  al.  2006a, b).  For  these  objects,  there¬ 
fore,  we  extended  the  Bergeron  grid  using  the  DO  models  from 
Hugelmeyer  et  al.  (2006)  as  a  guide.  However,  this  procedure 
likely  introduces  magnitude  errors  of  approximately  0.5  mag. 

Finally,  we  had  to  apply  a  reddening  correction  to  our  calcu¬ 
lated  absolute  and  observed  magnitudes  in  order  to  obtain  dis¬ 
tances  to  our  objects.  Since  we  don’t  know  how  much  of  the 
measured  interstellar  reddening  is  really  between  us  and  our  in¬ 
dividual  target  stars,  we  employed  the  same  iterative  reddening 
approach  as  described  in  Harris  et  al.  (2006),  assuming  zero  red¬ 
dening  at  100  pc  and  full  reddening  at  250  pc. 


describes  the  maximum  volume  in  which  a  star  can  contribute 
to  the  FF: 

An  t  'x  2  \  —  w 

Fmax  =  (rmax  —  rmin  j  e  ~5° 

where,/!  is  the  ratio  of  the  SDSS  DR4  spectral  observation  cov¬ 
erage  of  the  sky  (4783  degrees;  Adelman-McCarthy  et  al.  2006) 
to  the  total  sky  surface.  We  accommodate  the  hot_std  incom¬ 
pleteness  by  introducing  the  chs  factor.  Finally,  the  e~A  fac¬ 
tor  takes  into  account  the  exponentially  decreasing  WD  density 
with  distance  z  (in  parsecs)  from  the  Galactic  plane,  assuming  a 
Galactic  scale  height  of  250  pc  as  in  Harris  et  al.  (2006). 

The  chs  incompleteness  factor  is  not  a  constant  number  for 
all  WDs  in  our  sampled  due  to  the  additional  isolation  criteria 
placed  on  hot_stds.  The  hot_std  incompleteness  is  a  function  of 
both  Galactic  latitude,  b  and  object  brightness,  g.  At  lower  b, 
the  increased  number  of  stars  will  make  the  hot_std  unblended 
criteria  harder  to  meet,  thereby  increasing  the  incompleteness. 
Similarly,  brighter  objects  cover  more  pixels  on  the  CCD,  and 
thus  are  more  likely  to  be  thrown  out  of  the  hot_std  category 
due  to  the  blending  or  edge  criteria.  Eisenstein  et  al.  (2006a) 
estimates  that  hot_std  incompleteness  varies  from  0.55  to  0.60 
for  30  <  b  <  40  and  from  0.60  to  0.63  when  40  <  b  <  70,  then 
from  0.40  to  0.55  for  g  magnitudes  15  <  g  <  17,  and  from  0.55 
to  0.63  for  17  <  g  <  19.5  mag. 

Because  we  need  to  know  the  incompleteness  as  both  b 
and  g  vary,  we  queried  the  SDSS  database  to  obtain  a  full  two- 
dimensional  incompleteness  map.  We  formed  a  6  x  18  grid  with 
g  varying  from  14-19  in  1  mag  bins  and  b  varying  from  -90 
to  +90  in  10  degree  bins.  Our  method  was  to  determine  the  ra¬ 
tio  of  hot_stds  observed  spectroscopically  to  the  total  number  of 
potential  hot_stds  that  passed  the  photometric  criteria,  but  not 
the  isolation  requirements  for  each  of  our  ( g ,  b)  cells.  The  SDSS 
software  pipeline  separates  blended  objects  into  a  parent  (de- 
blended)  object  and  its  constituent  children.  We  allowed  the  par¬ 
ent  object  in  our  search,  but  not  the  children.  We  also  removed 
all  moving  objects  (e.g.,  asteroids)  and  objects  too  close  to  the 
edge  of  the  CCD  chips  to  make  sure  that  the  photometry  was 
properly  determined.  Table  1  shows  the  results  of  our  cus(9,t>) 
calculations.  Values  in  cells  are  given  in  [%]  and  denotes 
cells  with  no  spectroscopic  hot_stds  in  them.  Magnitude  bins 
are  1 .0  mag  wide  in  g  starting  from  given  magnitudes  in  the  first 
column.  Fatitude  bins  are  10  degrees  wide  extending  from  the 
values  in  the  first  row. 

Altogether  we  found  5778  photometric  (PH)  hot_standards 
and  1984  spectroscopically  (SP)  observed  objects  formally  clas¬ 
sified  by  the  pipeline  as  a  hot_std.  The  incompleteness  factors  of 
spectroscopic  hot_stds  for  each  cell  can  be  calculated  as: 


0  ri 

chs,  =  -jpjj  ■  footprint  -  100% 


where  footprint  (~  1 .39)  is  the  photometric  to  spectroscopic  cov¬ 
erage  ratio  of  the  sky  surface  in  SDSS  DR4.  A  detailed  map  of 
hot_std  incompleteness  for  SDSS  survey  is  provided  in  Table  1. 
The  dash  signs  in  the  table  denotes  cells  where  no  spectroscopic 
hot_stds  were  found,  which  does  not  necessarily  mean  that  there 
are  no  photometric  hot_stds  in  those  bins. 


3.2.  Completeness 

To  normalize  our  LF,  we  employed  the  standard  ( 1/Vmax ) 
method  as  introduced  by  Schmidt  (1998).  The  quantity  Vmax 


3.3.  Results 

Using  the  incompleteness  map  given  in  Table  1,  we  calcu¬ 
lated  our  LF  for  Mbo  1  bins  of  0.5  mag  width.  The  resulting 
luminosity  function  is  presented  in  Fig.  2  along  with  that  of 
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Table  1.  SDSS  DR4  spectroscopic  hot_std  sample  incompleteness  map. 
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Fig.  2.  Our  primary  WD  LFs.  The  solid  line  includes  all  stars  in  our 
sample  while  the  dotted  line  excludes  the  more  uncertain  magnetic  and 
binary  WDs.  The  two  dashed  lines  show  the  effect  of  different  tempera¬ 
ture  cuts:  our  preferred  cut  at  23  500  K  and  one  cut  further  at  20  000  K 
to  provide  more  overlap  with  the  Harris  et  al.  (2006)  sample,  shown  as 
the  solid  line  with  open  circles.  The  top  axis  shows  the  approximate  ef¬ 
fective  temperature  in  kilo- Kelvin  corresponding  to  plotted  x-axis,  Mb0i. 


Harris  et  al.  (2006).  It  covers  an  absolute  bolometric  magnitude 
range  between  -0.75  and  7.25  and  provides  some  overlap  with 
the  cooler  LF  of  Harris  et  al.  (2006). 

To  determine  exactly  which  bins  were  affected  by  our  tem¬ 
perature  cut,  we  calculated  the  LF  for  two  different  cutoff  tem¬ 
peratures,  our  nominal  23  500  K  and  one  at  20  000  K,  as  dis¬ 
cussed  in  Sect.  2.  The  two  LFs  converge  at  Mboi  ~  6.75,  which 
we  therefore  take  as  the  coolest  complete  bin  in  our  LF.  There 
appears  to  be  general  agreement  with  Harris  et  al.  (2006),  al¬ 
though  there  aren’t  any  complete  bins  in  common. 

Note  the  plateau  seen  in  our  LF  between  Mboi  ~  0.5-3.75 
(Tefr  ~  50000-100000  K).  Eisenstein  et  al.  (2006a)  point  out 
a  systematic  overestimation  of  Tes  in  their  model  fits  for  the 
hottest  WDs,  starting  for  DAs  around  7dr  ~  30000  K.  At 
T eff  ~  50  000  K,  approximately  A7hoi  =  4.0,  this  overestimate  is 
a  roughly  10%  effect.  Furthermore,  at  these  temperatures,  non- 
LTE  (NLTE)  effects  become  important  whereas  the  Eisenstein 
et  al.  (2006a)  fits  were  based  only  on  LTE  models.  To  correct  for 
these  two  mis-estimates,  we  applied  a  parabolic  temperature  cor¬ 
rection  determined  by  the  Eisenstein  et  al.  (2006a)  comparisons 
to  existing  literature  fits  and  NLTE  corrections  from  Napiwotzki 
et  al.  (1999).  The  resulting  LF  is  shown  in  Fig.  3.  As  expected, 
the  cool  end  of  this  corrected  LF  remains  largely  the  same,  but 
the  hot  end  changes  considerably,  the  parabolic  temperature  ad¬ 
justments  having  the  biggest  effect. 

With  these  corrections  included,  41  stars’  temperatures 
dropped  to  below  our  23  500  K  cutoff.  Table  2  lists  the  LF 
bin  values  both  before  and  after  these  corrections  (as  shown  in 


Fig.  3.  The  same  as  Fig.  2  but  with  reu  and  NLTE  corrections  applied. 


Table  2.  Number  of  all  types  of  WDs  per  magnitude  bin. 
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Figs.  2  and  3).  Columns  2  and  3  list  the  number  of  DAs  in  each 
bin  before  and  after  the  corrections.  Since  the  corrections  to  the 
DAO,  DB,  DO,  and  PG  1159  stars  were  not  needed,  they  are 
listed  separately  in  Cols.  4  through  7,  respectively.  Finally,  Col.  8 
lists  the  total  number  of  WDs  of  all  types  in  each  bin. 

As  Fig.  3  shows,  the  plateau  seen  in  Fig.  2  narrows,  but 
remains  even  after  the  NLTE  and  fit  temperature  corrections 
were  applied.  Figure  4  repeats  the  LFs  shown  in  Fig.  2,  but 
with  1 .0  mag  bins  to  increase  the  signal  in  each  bin;  again  the 
plateau  remains.  While  there  are  still  problems  with  the  hottest 
DA  WDs  temperature  determination  (see  conclusions),  a  glance 
at  Table  2  shows  that  the  DAs  start  appearing  at  Mboi  ~  2.0, 
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Fig.  4.  Our  WD  LF  as  in  Fig.  3,  but  with  1.0  mag  bins. 


Fig.  5.  DA  (solid  circles)  and  non-DA  (open  squares)  LFs  plotted  with 
the  PG  survey  DA  LF  (open  circles)  of  Liebert  et  al.  (2005).  The  mag¬ 
nitude  bin  width  is  1 .0  and  the  bins  are  centered  on  integer  magnitudes. 
Thinner  line  right  below  SDSS  DA  LF  is  the  DA  LF  without  magnetic 
DAs  and  binary  DA+M  systems.  In  case  of  non-DA  LF  removing  mag¬ 
netic  and  binary  DB+M  systems  from  the  WD  sample  affects  the  LF 
near  Mb0i  =  7  and  the  rest  of  the  thin  line  overlaps  with  the  thick  one. 
The  top  axis  shows  the  approximate  effective  temperature  in  kilo-Kelvin 
corresponding  to  plotted  x-axis, 


precisely  where  our  plateau  begins.  It  appears  we  are  seeing  the 
emergence  of  the  non-DAs  into  the  LF. 

4.  DA  and  non-DA  luminosity  functions 

Table  2  shows  that  the  hot  end  of  our  LF  is  made  up  entirely  of 
non-DA  WDs.  We  therefore  divided  our  main  LF  into  two  LFs, 
one  with  only  DAs  and  one  for  non-DAs.  We  kept  the  23  500  K 
TeS  cutoff  for  the  DA-only  LF,  but  were  able  to  decrease  the 
temperature  cutoff  to  16  000  K  for  the  non-DA  LF.  Using  bins 
of  1 .0  mag  centered  on  half  integer  and  integer  magnitudes,  we 
plot  the  resulting  new  LFs  in  Figs.  5  and  6,  respectively.  We  also 
compare  to  the  DA  LF  of  Liebert  et  al.  (2005),  after  converting 
it  from  My  into  Mbd  space  densities. 

The  first  thing  to  notice  is  that  our  new  DA  LF  and  that  of 
Liebert  et  al.  (2005)  are  in  very  good  agreement.  The  hot  end 
of  the  non-DA  LF  has  an  interesting  shape,  with  a  sudden  drop 
in  the  LF  near  ;Uhoi  =  2.  A7hoi  =  2,  however,  this  is  precisely 
where  the  DAs  start  appearing.  If  we  follow  the  initial  slope  of 


T,„  [kK] 


Fig.  6.  The  same  as  in  Fig.  5  but  bins  are  centered  on  half  integer  mag¬ 
nitudes.  In  both  plots,  we  can  see  the  appearance  of  the  first  DAs  affect 
a  dip  in  the  non-DA  LF. 

Table  3.  DA  and  non-DA  WD  distribution  per  magnitude  bin. 


Teff 

[kK] 
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magnitude 
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9 
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6.0 
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21 
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7.0 

198* 

18 

11 

24 

8.0 

- 

16 

19 

9.0 

- 

72 

15 

10.0 

114* 

the  non-DA  LF  and  extend  it  before  the  drop-off  to  Mboi  =  2,  it 
leads  directly  to  the  start  of  the  DA  LF.  Thus,  it  seems  this  drop 
off  is  simply  due  to  the  appearance  of  DAs.  These  luminosity 
functions  represent  the  first  time  this  separation  of  the  WD  pop¬ 
ulation  has  been  seen.  Table  3  provides  the  numbers  behind  the 
separate  DA  and  non-DA  LFs  as  shown  in  Fig.  6.  The  last  DA 
and  non-DA  bins  (marked  with  asterisks)  are  incomplete  due  to 
our  temperature  cutoffs  at  23  500  K  and  16  000  K  respectively. 

5.  Conclusions 

As  shown,  our  WD  LF  agrees  very  well  with  previous  efforts, 
but  extends  to  significantly  hotter  stars  and  shows  clear  emer¬ 
gence  of  DAs  and  non-DAs  in  the  general  WD  cooling  curve. 
However,  due  to  the  problems  with  hot  DA  WD  temperature  de¬ 
termination  the  DA  and  non-DA  space  densities  cannot  be  com¬ 
pared  directly  at  the  hot  end  of  these  LFs.  Although  the  sys¬ 
tematic  DA  Te ff  overestimation  of  Eisenstein  et  al.  (2006a)  WD 
catalog  data  has  been  corrected  and  the  pure  hydrogen  NLTE 
vectors  have  been  applied,  yet  there  are  more  severe  problems 
affecting  the  hottest  DA  WD  temperature  and  log  g  determina¬ 
tion.  These  include  lower  spectra  quality  of  the  fainter  HSs  or  the 
presence  of  metals  in  the  hot  DA  WD  atmosphere  visible  as  the 
onset  of  the  Balmer  line  problem  (Vennes  1999)  in  the  spectra. 
Barstow  et  al.  (2003)  and  Vennes  et  al.  (2005)  investigated  the  ef¬ 
fect  of  metallicity  influence  on  the  hot  DA  WD  (above  50  000  K) 
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temperature  determination.  They  have  shown  that,  due  to  heavy 
elements  present  in  the  DA  atmosphere,  in  some  cases  the  hot 
DA  temperatures  can  be  overestimated  by  at  least  20%,  which 
translates  to  ~1  Mboi  more  luminous  WDs.  In  such  circum¬ 
stances  the  space  density  of  the  DA  WDs  seems  reasonably  well 
established  for  M\,a\  >  3  mag  and  is  not  certain  for  the  more  lu¬ 
minous  part  of  the  DA  LF.  In  that  regard,  the  non-DA  LF  is  likely 
the  most  reliable  out  of  both  our  LFs.  It  covers  -0.5  to  9.5  bolo- 
metric  magnitude  range  and  shows  the  drop  off  near  ;Whoi  =  2 
which  we  interpret  as  the  non-DA  to  DA  atmosphere  conversion 
(Figs.  5,  6).  We  can  see  that  in  the  DB  gap  temperature  range, 
30000-43000  K  5-6),  our  non-DA  LF  has  only  a  slight 
space  density  decrease,  suggesting  the  number  of  recently  dis¬ 
covered  hot  DB  stars  (Eisenstein  et  al.  2006a, b)  has  nearly,  if 
not  completely,  filled  the  “gap”. 

Note  that  the  error  bars  presented  in  the  figures  are  formal 
errors,  calculated  as  a  sum  of  l/(Vniax)2  for  each  luminosity  bin, 
and  they  do  not  take  observed  or  absolute  bolometric  magnitude 
error  determinations  into  account.  In  fact  the  Fmax  method  un¬ 
derestimates  LF  error  bars  by  a  factor  of  10  (Geijo  et  al.  2006) 
near  the  9-7  Mboi  magnitude  range  where  our  LF  has  its  cool 
end,  and  in  that  regard  the  agreement  between  our  LF  and  others 
can  have  larger  tolerance.  In  their  theoretical  work  Geijo  et  al. 
(2006)  did  not  make  any  LF  error  estimations  for  the  hotter  re¬ 
gions  our  LF  covers,  so  we  cannot  determine  the  error  quality 
in  that  range,  but  we  can  expect  the  ones  we  present  are  also 
underestimated. 

In  short,  we  have  calculated  a  new  WD  LF  with  an  ex¬ 
tended  range  of  luminosities  based  on  the  well-understood  and 
explored  SDSS  WD  data  set.  However,  to  improve  the  quality 
of  the  LF  at  its  hot  end,  new  follow  up  observations  of  low  sig¬ 
nal  to  noise  spectra  are  necessary  as  well  as  careful  analysis  of 
the  metal  content  in  the  WD  spectra.  Further  exploration  of  our 
LF  should  contribute  a  number  of  new  insights  into  early  white 
dwarf  evolution. 
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